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Abstract 

Elastic properties of soft, three-dimensional dimers, interacting through site-site 
n-inverse-power potential, are determined by computer simulations at zero tem¬ 
perature. The degenerate crystal of dimers exhibiting (Gaussian) size distribution 
of atomic diameters - i.e. size polydispersity - is studied at the molecular number 
density l/\/2; the distance between centers of atoms forming dimers is considered 
as a length unit. It is shown that, at the fixed number density of the dimers, in¬ 
creasing polydispersity causes, typically, an increase of pressure, elastic constants 
and Poisson’s ratio; the latter is positive in most direction. A direction is found, 
however, in which the size polydispersity causes substantial decrease of Poisson’s 
ratio, down to negative values for large n. Thus, the system is partially auxetic for 
large polydispersity and large n. 
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1 Introduction 


Among the characteristics of real materials, elastic properties are the ones that 
are crucial from the point of view of various applications. Recently, increasing 
interest is observed in the area of materials with unusual elastic properties - 
an example are materials with negative Poisson’s ratio . When stretched 
(or resp. compressed) such materials increase (resp. decrease) their size not 
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only in the direction of the stretch (resp. compression) but also in one or 
more directions transverse to it [^. Materials exhibiting such counterintuitive 
behavior are known as anti-rubber [^, dilational materials or auxetics j^. 
In this paper, the last name will be used for its brevity. 


As for isotropic systems the Poisson’s ratio is direction independent, any 
isotropic material can be either auxetic or not. For anisotropic systems, which 
are the subject of the present paper, the situation is more complex. Apart 
from the alternative that the Poisson’s ratio is in all directions negative or 
not, a third possibility exists: in some directions it is negative and in other - it 



and nonauxetics correspond, respectively, to negative or nonnegative Poisson’s 
ratio in all directions. Materials with Poisson’s ratios negative in some direc¬ 
tions and nonnegative in the other ones will be referred to as partially auxetic 
materials or just partial auxetics. 


The very existence of auxetics and partial auxetics encourages one to study 
influence of various micro-level mechanisms on the elastic properties of con¬ 
densed matter 13l444l|. An important problem in this context is the influ¬ 


ence of disorder on the Poisson’s ratio. For some two-dimensional (2D) mod¬ 
els which are elastically isotropic (at least for small deformations) and for 
some anisotropic 2D models, it has been found that various kinds of disorder 
(such as vacancies, interstitials, aperiodicity or polydispersity) introduced in 
crystalline structures, increase Poisson’s ratio 45l-l50|. Studies of some three- 
dimensional (3D) molecular models revealed similar effect caused by va¬ 
cancies, what might suggest a general conclusion that (at least in simple molec¬ 
ular systems) disorder weakens or eliminates auxeticity. 


Recently, however, it has been found that in one of the simplest 3D systems - 
soft spheres, forming the f.c.c. phase, there exists a direction in which one of 
the form of disorder (polydispersity) decreases Poisson’s ratio, which reaches 
negative values down to -1 |52|. This observation might hnd practical appli¬ 
cations, e.g. in production of partial auxetics of latex meso- or macroscopic 
spheres for which size polydispersity is easy to control. Unfortunately, the 
mentioned lowest Poisson’s ratios were observed only for the softest spheres 
with n = 6, i.e. when the f.c.c. structure may be unstable with respect to the 
b.c.c. structure 53 . 


Taking into account that elastic properties of many-body systems strongly de¬ 
pend on details of the intermolecular interaction potential and, in particular, 
on the molecular shape, it is interesting to study other 3D systems searching 
for stable structures for which increasing polydispersity can lead to more neg¬ 
ative Poisson’s ratios in some directions. As it is meaningful to start with the 
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simplest molecular shapes, the natural candidate for a hrst 3D molecule to 
investigate is a polydisperse dimer. 


The aim of this study is to determine the influence of disorder introduced 
to a solid of soft dimers by the size dispersion of their atoms, on the elastic 
properties of the model. The study concerns the zero temperature, T = 0, and 
the number density N/V = 1/\/2, where V is the volume and N is the number 
of dimers. At this density atoms of dimers can be arranged in a perfect f.c.c. 
lattice of lattice constant equal to the distance of atomic centers in the dimer. 
The dimers form a degenerate crystalline (DC) phase - an analogue of that 
observed in 2D The phase is characterized by non-periodic positions and 
orientations of molecules, whereas the atoms forming molecules are arranged in 
a perfect f.c.c. lattice (when the atoms are identical) or (for slightly different 
sizes of atoms) in a lattice close to the f.c.c. one. The DC phase is known 
to be thermodynamically stable and effectively cubic from the point of view 
of elastic constants for hard dimers and is expected to preserve those 
properties for soft dimers with short ranged interactions which are studied in 
the present paper. 


The soft dimer system can be seen as a rough model of some man-made 
nano-, meso- or macro-particle systems, e.g. systems of latex particles j^. 
Recently, an increasing interest is observed in systems of soft particles 


a. 


especially in the context of physics of colloids. This is due to the fact that 
the latter have a broad range of possible applications. Systems of soft poly¬ 
disperse particles seem to be more realistic models of colloids than those built 
of identical particles. Polydispersity is also important for nanoparticle crys¬ 
tals (m-^, where the speed of crystallization and resulting crystal quality is 
strongly related to size dispersion of particles. 


Studies of elastic properties of models of polydisperse systems can be thought 
of as a hrst step on a path to better understanding mechanisms and phenomena 
that govern the elastic properties of real systems. Such a knowledge is not only 
important from the point of view of fundamental research but also useful from 
the material engineering point of view - it can help to engineer materials with 
desired elastic properties. 


2 Description of the studied model 


The model studied in this work consists of N diatomic molecules (dimers) in 
periodic boundary conditions. The dimers are rigid i.e. the distance between 
centers of the atoms forming each one is constant and equal to a (which 
constitutes a unit of length) for each dimer. Initially the atoms form the perfect 
f.c.c. lattice and the molecules form a degenerate crystalline (DC) phase. The 
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atoms interact via purely repulsive potential of the form: 



n 


^ij ( dij Uq 


( 1 ) 


where dij = {di + dj)/2, di, dj are the diameters of the interacting spheres 
(atoms), Vij is the distance between their centers and Uq is the energy unit. 
The exponent n is further referred to as the hardness parameter and its in¬ 
verse, 1/n, as the softness parameter. (These names come from the observation 
that when n ^ oo the interaction potential described by Eq. ([I]) tends to the 
hard-body potential, which is widely used in computer simulations and con¬ 
densed matter theory j^.) Atomic interactions are assumed to be short-range 
ones, i.e. only the nearest neighbors (shearing a face of their Voronoi polyhe- 
dra) interact. Since the molecules are considered as rigid, interactions between 
atoms forming a molecule are neglected. 

2.1 Atomic size polydispersity 

In the studied system, atomic diameters were generated randomly according 
to the Gauss distribution function with a given standard deviation 6 



( 2 ) 


In above equation, d is the atomic diameter and a = (d). Further, in this 
work, 6 is treated as the measure of disorder in the system and is referred to 
as the polydispersity parameter. Generated values of atomic diameters were 



Fig. 1. Example of studied structures with atomic size polydispersity. The structure 


shown consists of 432 molecules and has the polydispersity parameter 5 equal to 3%. 
In (a) the atoms are drawn in the shades of gray, indicating their size relative to 
other atoms. The ones drawn in white are the smallest and those drawn in black are 
the largest. The gray cube indicates the periodic box. In (b) connections between 
atoms in the same structure as in (a) are shown. 
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randomly distributed among atoms in the structure. In the Fig. [T] examples of 
studied structures are presented. It can be seen that even in the presence of 
atomic size polydispersity atoms form nearly perfect f.c.c. lattice. 


3 Basic formulae 


In the following, it will be assumed that the studied system of soft dimers 
exhibits effective regular symmetry. For the latter symmetry, when the exter¬ 
nal pressure p is zero, elastic properties can be described by just three elastic 


const ants •, Cxxyy and Cxyxy 

nient [ 


62| . However, when p 7 ^ 0, it is more conve- 


45| to use other elastic constants, Bijki, obtained by expanding the 


free enthalpy (Gibbs free energy), G, with respect to the components of the 
Lagrange strain tensor, e: 


AG/Kef = A(F 


pV) ^ ( BjjkiEjjEki . 

ijkl 


(3) 


F in the above equation is the (Helmholtz) free energy and G-ef is the volume of 
the deformation-free (reference) system. The elastic constants Bijki are simply 
related with the constants Cijki |^, 0 


Bijki Cijki P i^^ik^jl T ^il^jk ^ij^kl) 


(4) 


where 6ij is the Kronecker delta, equal to 1 for i = j and zero otherwise. 

In the case of cubic symmetry, considered in this work, only the following 
elastic constants are of interest (the Voigt notation is used below (^): 


Bn = Cn — P , 

(5) 

Bi2 = Ci2 + p , 

(6) 

= Cee — P , 

(7) 


and the free enthalpy expansion (|3]) takes the form: 

AG/Vef = -Bii (^elx + e^y + 

B\2 i^^xx^yy ^yy^zz “ 1 “ ^zz^xx) (^^xy ^yz ^zx'^ * 

For cubic systems, instead of constants defined by Eq. 0 - 0 . the below elastic 
moduli are also often used: 
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(9) 

( 10 ) 

( 11 ) 


= -Bee 5 
1 / 

^'■2 = 2 (-Bn — -B 12 ) , 

where B is the hulk modulus, and the /ij are the shear moduli. In the following, 
both sets of the quantities will be used while presenting and discussing the 
results. 

The Poisson’s ratio is dehned as the negative ratio of transverse to longitudinal 
strain when only the longitudinal component of the stress tensor is inhnites- 
imally changed. For anisotropic systems it depends, in general, on both the 
longitudinal and the transverse direction. For cubic symmetry, the Poisson’s 
ratio along two high-symmetry longitudinal directions ([100] and [111]) does 
not depend on the choice of the transverse direction [l^ and can be expressed 
by the elastic constants, Bij, in the form j^: 


Bi2 1 SB — 2/i2 
^ Bu + 5i 2 ^ 2 35 + /i2 ’ 

Bii 2 Bi2 — 25(56 1 35 — 2/ii 

^ ^ 2 (5ii -|- 25 i2 + Bqq) 2 SB ni 


( 12 ) 

(13) 


In this paper, the Poisson’s ratio along the longitudinal direction [110] and 
two transverse directions [110] and [001] was also computed: 


^[110] [ilo] 


B\i + 5ii5i2 — 25^2 ~ 25ii566 
Bii -|- 5 ii5i2 — 25^2 + 25ii566 
35(3/i2 - /ii) - 4/iip2 
35 (3yU.2 + /ii) + Afiifi2 


(14) 


_ 45i2566 

^[110][001] j2)2 \ tz) jd 9 r2 i 9 R R 

-Dll + -D1I-D12 — ^-Di 2 + ^-Dll-Dee 

65yUl — 4/ilyU,2 

35 (3yU.2 + /il) + 4yUl/i2 


(15) 


4 Numerical determination of the elastic properties 


In order to obtain the elastic constants, a sample was considered which, at 
equilibrium, is described by three vectors uq, 6o and Cq parallel to the unit 
cell edges of the f.c.c. lattice. Components of those vectors form columns of a 
reference box matrix (o^ . 67 ] H. Applying a uniform deformation that trans¬ 
forms the sample into a parallelepiped described by the box matrix h one can 
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( 16 ) 


write the Lagrange strain tensor as j66l . l67 


£ = (h'”^ h' h - l) /2 , 

where h' is the transposed matrix h and I is the unit matrix. Differentiating 
the free energy F with respect to the strain tensor components, one obtains 
all the elastic constants and pressure. Is is worth to note, that in the case of 
zero temperature (T = 0°K), the free energy is just equal to the potential 
energy of the system. 


After determining equilibrium positions and orientations of all the molecules 
in the reference state for a given 5, the following deformations were applied: 

(Iff ^a) HoipSoip , 

for which new equilibrium positions and orientations were determined. In fll7l) 
a, f3 indicate directions, and is a small real number. The relations below 
define the elastic constants and pressure: 


dF 


d‘^F 


dil 






d^F 


diadip 


ia=i/3=0, a^P 


= -p, 

= -Bii , 

= - P ■ 


(18) 

(19) 

( 20 ) 


The ilge constant was obtained by applying the deformation of the form: 


'^ 0 ( 0 : -^-^ 0 : 0 : ? 

^ap Hqi^p T CapHaa i (21) 

for which new equilibrium positions and orientations were determined. In fl21l) 
Cap = Cpa is a small real number. For Fl^x = Hyy = H^z this leads to: 

= 6(2rB66ffp) • (22) 

Cc.^=o 


02 F 

Wp 


5 Results 


In the Figs. [2] (a) and [2] (b) values of the pressure and the bulk modulus 
are shown, respectively. Is is clearly visible that increasing polydispersity and 
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Fig. 2. The values of (a) pressure and (b) bulk modulus B for polydisperse dimers, 
plotted against the softness parameter (1/n). The thick solid line represents values 
for the system with (5 = 0. 


increasing n leads to increase of those quantities as it was earlier observed in 
other systems 

The observed behavior can be qualitatively understood by noticing that dense 
packing of identical molecules in the hard body limit requires less space than 


dense packing of polydisperse molecules of the same average size j^. As all 


the considered systems are studied at the same volume, their densities relative 
to close packing grow with increasing polydispersity. Increase of the relative 
density implies an increase of the pressure and the bulk modulus. Obviously, 
the latter quantities show the largest increase for the largest n, when the 
interaction potential is the steepest one. 


In Fig. [3] the elastic constants Bij are shown. It can be seen that all elastic 
constants also increase with increasing polydispersity and increasing vaule of 
n. Presented behavior of elastic constants is also similar to that observed for 
other 2D and 3D systems 47H50l. l52| and can be qualitatively explained by 


the close packing argument used above. 





Fig. 3. The values of elastic constants (a) Bn, (b) Bn and (c) Bqq. plotted against 
the softness parameter (1/n). The thick solid line represents values for the system 
with (5 = 0. 


















In Fig. lDplots oi Hi!B ratios are presented. As in the case of earlier studies j47- 
50 . 5^ for both 2D and 3D, this ratios seem to tend to 0 for any nonzero 
polydispersity, in the hard interactions limit {n —>■ cxo). 




Fig. 4. The values of h/B ratio for (a) yi and (b) H 2 -, plotted against the softness 
parameter (1/n). Thick solid line represents values measured for the system with 


<5 = 0. 


According to the formulae (IT^ and (IT^ the mentioned above asymptotics of 
Hi/B implies that the system behaves like rubber {^rubber = 0.5) when a small 
stress is applied along the directions [100] and [111] in the limit n —)> oo. This 
is shown in Fig. El where the plots of Poisson’s ratio are presented. It can be 
seen there that, typically, presence of any nonzero polydispersity causes an 
increase of the Poisson’s ratio in most of directions corresponding to the high 
symmetry axes. It is worth to notice that in the cubic system studied, the 
Poisson’s ratio i^[iio][ooi] exceeds the maximum value 1/2 allowed for isotropic 
materials. One should also stress that a pair of exceptional directions (see 
Fig. [Sl[d)) exists in the system for which the Poisson’s ratio decreases with 
increasing polydispersity. This occurs when the direction of the deformation 
is [110] and the response of the system is measured in the direction [110]. For 
this pair of directions, the Poisson’s ratio reaches negative values both for small 
and large n limits! The decrease for large n is, however, much more signihcant. 
This new behavior of Poisson’s ratio is opposite to what was discovered in 
the soft polydisperse spheres system j^, where the most negative values of 
Poisson’s ratio were observed for small n only. 


The observed dependence of surprising in the context of results 

obtained for various 2D isotropic systems including the DC phase of 2D 
dimers is of interest from the point of view of production of partial aux- 


etics. As shown in Fig. |5](d), solid structure of dimers obtained by ‘gluing’ 
together pairs of polydisperse soft spheres arranged in the f.c.c. lattice will 
show the Poisson’s ratio as low as = —0.8. 
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b) 


1/n 


V« d) 


Fig. 5. Poisson’s ratio plotted against the softness parameter (1/n). Plots 
(a), (b) show the Poisson’s ratio measured perpendicularly to the direction of the 
high-symmetry axis (a) [100] and (b) [111]- Plots (c), (d) show Poisson’s ratio per¬ 
pendicularly to the direction [110] in the direction (c) [001] and (d) [110], respec¬ 
tively. 


6 Conclusions 


The elastic properties of polydisperse soft dimer system in three dimensions 
were studied by computer simulations in the static (zero temperature) case. 
Disorder was introduced to the DC f.c.c. structure in the form of atomic size 
dispersion with polydispersity parameter 6. The influence of polydispersity on 
elastic properties of 3D DC dimers system was determined. It was found that 
increasing polydispersity and hardness of the interaction potential causes a 
signihcant increase of elastic constants and pressure. 

The simulations have shown that Poisson’s ratio of the dimer system measured 
in [100] and [111] directions increases to its maximum positive value, 1/2, in 
the presence of any nonzero polydispersity and sufficiently hard interaction 
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potential. An increase of Poisson’s ratio was observed also in other directions. 
Such a behavior was earlier observed in the case of 2D 
systems and seems to be a typical result. 


47-50 and 3D 52 


It is worth to stress, however, that for the longitudinal direction [110] and 
transverse direction [110] the Poisson’s ratio decreases reaching negative Ydlwes 
along with increasing polydispersity value both for large and small n. This 
unusual phenomenon is of interest from the point of view of manufacturing 
stable auxetics of cubic symmetry. 


The questions whether the effect of auxeticity enhancement by polydispersity 
is characteristic for the f.c.c. lattice only or if it can be observed for elastic 
symmetries other than the f.c.c. one, e.g. for packing of spheres in the h.c.p. 
lattice or stacking fault lattices, will be a subject of a separate paper. Studies 
of elastic properties of other molecular shapes are in progres. 
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